Vacuum fluctuations for spherical gravitational impulsive waves 
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Abstract We propose a method for calculating vacuum fluctuations on the background of 
a spherical impulsive gravitational wave which results in a finite expression for the vacuum 
expectation value of the stress-energy tensor. The method is based on first including a 
cosmological constant as an auxiliary constant. We show that the result for the vacuum 
expectation value of the stress-energy tensor in second-order perturbation theory is finite 
if both the cosmological constant and the infrared parameter tend to zero at the same rate. 



1. Introduction 

It has been known for a long time that plane waves do not give rise to vacuum 
fluctuations [1,2]. One may check whether the same result is true for spherical waves as 
well. One possible check is to use the Nutku-Penrose metric [3] for spherical impulsive 
gravitational waves and to study whether vacuum fluctuations arise when a scalar particle 
is coupled to the gravitational field through this metric. 

We applied this procedure to different holomorphic warp functions [4-8] and found no 
finite fluctuations in the Minkowski case. We must stress the fact that in these calculations 
first-order perturbation theory was used. This method gave a null result for the Minkowski 
case. Taking only first order perturbation theory was, perhaps , the weak point in these 
calculations. One conjecture [9] was that the null result would change in second-order 
calculations. Here our results on this problem , performed in second-order perturbation 
theory, will be reported. 

It is common knowledge that massless particle production is prohibited by the presence 
of conformal symmetry [10]. One may argue that this fact is not relevant here since the 
Nutku-Penrose metric is not conformally related to Minkowski space. This would be true 
if an exact calculation were made. One must note at this point that we are using only 
approximate methods to solve this problem. We are perturbing around the Minkowski 
space to calculate the fluctuations in the new space. This perturbation does not change 
the essential properties of the solutions of the Minkowski space, but just modulates them by 
multiplying by a factor. The Minkowski Green function is just a monomial over distance. 
We could not find a way to regularize it and still retain non-zero terms, since there are none 
left when the worst divergence is subtracted. Our perturbative calculation just modifies 
the Minkowski expression by multiplying it by a finite term. As a result we could not get 
any finite result in this case in first-order perturbation theory for the cases studied. 

If we had a mass parameter in the problem, we could get an expansion whose second 
or third term could give finite fluctuations. This was not possible here. We still hoped 
that the second-order perturbation calculation might give a signal for a possible way out 
of this impasse. 

One way to introduce a mass parameter without changing the Nutku-Penrose solution 
too much is to study a similar case in de Sitter space where one obtains an impulsive wave 
solution by a simple manoeuvre [11]. The scalar curvature of the de Sitter space has 
dimensions of mass squared. Here the calculation in first order gives finite results for all 
of the different choices of the warp functions studied [8,12]. If we carry the calculations to 
second order, at least for the special choice studied here, we get infrared divergences. We 
conjecture that this behaviour is generic for a wider class of warp functions, if not for all 
possible forms, since in first-order calculations we get the same general behaviour for all 
the different choices we have studied. We verify this conjecture for another warp function 
which has shown to be representative for a wider class of warp functions. 



We take these divergences for the signal we are looking for. The Green function we 
calculate includes terms which are inversely proportiona to the square of an infrared mass 
which should bectaken to zero at the end of the calculation. One natural way to get rid of 
these two parameters (the curvature of the de Sitter space and the infrared parameters) is to 
set them proportional to each other. In one sense both terms are introduced for technical 
reasons. We went to de Sitter space to perturb around a space which is different from 
the Minkowski in the first place, but one that still retains the Nutku-Penrose solution. 
The infrared parameter was introduced just to be able to calculate the Green function 
unambigiously. There were no massive parameters in the model initially, and they should 
not be present in the end result. Taking these two terms proportional to each other and 
then sending them to zero cancels these two auxiliary parameters completely, still retaining 
a finite contribution to < T^ v >. 

At the end of the calculation all reference to the de Sitter universe is gone and our 
result is what one would expect to get in the Minkowski case. We have the vacuum 
expectation value of one component of the stress-energy tensor, 
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where S(v) is the Dirac delta function. In the Nutku-Penrose metric the only non-zero 
component of the curvature tensor is proportional to the same function. One may an- 
ticipate, following Deser [1], that any non-zero result may be proportional to the same 
form. 

Our conclusion is that to obtain a non-zero result for the vacuum fluctuations in the 
case studied here, first we have to work in de Sitter space, and thus explicitly break the 
conformal symmetry in the metric we perturb about. When this calculation is carried to 
second order, infrared divergences are traded for the for the curvature of the space to get 
a finite result. At the end we are back in the metric which describes an impulsive wave in 
Minkowski space. 

In section 2 we describe the second-order calculation in Minkowski space. We pass to 
the de Sitter case just by multiplying the Minkowski result by a conformal factor. Before we 
go to the coincidence limit we set the two auxiliary parameters proportional to each other. 
We then take the coincidence limit and perform differentiations on the Green function to 
obtain a finite value for the vacuum expectation value of the stress-energy tensor. We 
conclude with a few remarks. 



2. Minkowski calculation 

We start with the Nutku-Penrose metric [3] 

ds 2 = 2dudv - u 2 \d( + v<d(v0f(()d(\ 2 . (1) 

Here v is the retarded time, u is similar to the radial distance and ( is the angle of the 
stereographic projection, f(() is the Schwarzian derivative of h(() which is the holomorphic 



warp fuction describing the impulsive wave, and Q(v) is the Heavyside unit step function. 
For different choices of the warp function, the non-zero component of the curvature tensor 
is multiplied by a different function, but the essential characteristics of the metric do not 
change. The Nutku-Pentrose solution corresponds to the snapping of a cosmic string, 
giving rise to a spherical impusive wave. 

In the past we calculated vacuum fluctuations by computing the stress-energy tensor 
of a scalar field in the background of this metric. Our warp functions h(() included a 
parameter which is related to the string parameter G mu rs 10~ 6 that we wrote as 5 or 

e. We used h(() = (^r§) for different values of (A,B,C,D) [4-8]. We got the null 

result in all of these cases if we perturbed around the Minkowski space. If we multiply 
this metric by a conformal factor, , (l + Ap) we get a spherical wave in de Sitter space 
[11]. Then there are finite fluctuations even in first-order perturbation theory that are 
proportional to the square of the scalar curvature [12]. 

To investigate the same phenomena in second-order perturbation theory, here we use 
a different choice of the warp function. We take h = e a< > . There are two reasons why this 
function is chosen. First, this is an important special case, first discussed by Gleiser and 
Pullin [13]. It has the special property that f(() is independent of zeta although /i(C), the 
decisive term in the solution, depends on the same parameter. This fact allows us to write 
our solutions as sums of functions of £ and (. Our essential motive for studying this case 
is not this technical point, though. A similar problem was solved exactly [14], in another 
context. There were somewhat vague indications of particle production. Here we want to 
see whether the vacuum expectation value of one component of the stress-energy tensor 
is proportional to the non-zero component of hte curvature tensor. Although the result 
concerning particle production would not allow us to conclude anything definite about our 
new calculation, we may still argue that an unambigious result should exist for our latter 
calculation as well. We will just use the exact result as a guard against possible infrared 
problems. If such problems arise in our expansion we will know that they are due to 
technical factors, since the exact solution does not have them, and we will try to regularize 
them. 

We start with h = e a<5 . This gives / = — %- resulting in a metric 

ds 2 = 2dudv - - (dx 2 (2u - va 2 Q(v)) 2 + dy 2 {2 + va 2 Q(v)) 2 ) . (2) 



Here ( = x + iy. If we write the d'Alembertian operator in this metric , we get 



for the exact operator. We multiply the d'Alembertian operator given in equation (3) by 
u 2 and expand the operator up to second order in a 2 : 



L 11 = 2u 2 d u d v + 2ud v -d 2 x -d 2 -^(d 2 x -d 2 )-a* {^-d u - V \ + l^(d 2 x + d 2 )^j . (4) 

We can construct the vacuum expectation value of the stress-energy tensor from the two- 
point function through differentiation, after the coincidence limit is taken and the ap- 
propriate regularization is done. The two-point function Gf is found by summing the 
eigenfunctions of the related Sturm-Liouville problem. 

In zeroth order we get the empty space solution from the solution </>(°) [5] : 

G (o )=V £|^V) = A 

F Y A (u-u')(v-v')-^f((x-x') 2 + (y-y') 2 )' {) 

Here A is a constant. The first order solution is written as 0^1) = 0^0)/. Here is the 
zeroth-order solution. The first order solution is the product of the zeroth-order solution 
and another function. This ansatz for 0^ is dictated by the differential equations, and is 
not just an ad hoc guess. We find that / just modulates the zeroth-order solution, and 
does not essentially change it. The signature of <j>(°> is seen in the ultra-violet behaviour 
ofGi? to a large extent. 

We find that / obeys the differential equation 



L 2 f=-(k 2 -k 2 ) (6) 
u 



where L 2 is defined as 



L - (-2iR— - 2i (k — k— \-— 92 2 92 l JLJL\ ( 7 ) 
\ ds \ 1 dx 2 dy J dx 2 dy 2 dsdv R dv J 

Here k\,k 2 ,R , and K are parameters of the zeroth-order solution, s = -. They are 
integrated over to get the two-point function Gq- 

To calculate / we make the ansatz / = vfi(s,x,y) + f 2 (s,x,y) as suggested by the 
form of equation (6). This ansatz yields two equations, partly coupled: 

L 3 ii=^ 2 2 -fc 2 ), (8) 

and 

T,n f n = ( 1. 

ds R 

where 



Uh = 64- + % ) h (9) 



These equations are simply integrated over by multiplying the right-hand side by the 
inverse of the operator L 3 . 

We end up with 
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To get Gf, we have to calculate 0[f] where the operator O is given by 
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When we perform this operation, we get Gf equal to 

[(./• - *') 2 - (y - y) ] \ A, \ s _ s , ){ ] 

9(,) S 2 + 9K)/ » 3 9(t;)-^9(t/) . 
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where [ ] = (u — u'){v — v') — ((x — x') 2 + (y — y') 2 ), and Ai,A 2 ,A$ are constants. 
We see that this result is of the Hadamard form. No problems seem to arise in the infrared 
region. We find that all these terms have the same type of ultraviolet singularity as the 
flat part. We could not find a finite part of this expression. 

If we go one order higher, we end up with the differential equation 
L 2 g = v 2 (^ + 3s 2 (k 2 + k 2 ) + (k 2 - k 2 ) 2 

f-s iKs 2 (k 2 -k 2 ) 2 s s iK(K 2 -k 2 ) 2 s 4 \ 
+V {^ + ^R- AR 2 + 2A& ) < 14 > 

when we make the ansatz 4> 2 = <Po9- We take g = v 2 gi(x, y, s) + vg 2 {x, y, s) + gs (x, y, s). 
Going through similar steps as described above, we find 

s 2 is^ki + kp (fcj-fcpy 

9l ~ — + 2R + 32R 2 ■ (15) 

All these terms give two-point functions, Gp, in the Hadamard form. All are finite in the 
infrared sense. 



We use gi to get g 2 : 
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At this point we start getting problems. All the terms in this expression are divergent in 
the infrared region. When we apply the operator ) to g 2 and take x = x' , y = y', we get a 
term proportional to 

-i(u — u')(v — v') 
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We can start with a massive scalar particle and set the mass equal to zero at the end of 
the calculation. Then the above integral reduces to H (m^/ (u — u')(v — v')) where H is 
the Hankel function of order zero. It degenerates to a logarithm when m tends to zero. We 
have to take derivatives of Gi? to find the vacuum expectation value of the stress-energy 
tensor. The term with m decouples if we differentiate the logarithm with respect to u or 
v. The finite part of < > will not have logm if there is a finite part. We call this type 
of divergence 'harmless', in this sense. 



For# 3 , we get 
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We see that when x = x', y = y' all these terms give rise to expressions that are 
proportional to 
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which are linearly divergent at the upper limit . If we use a massive field as an infrared 
cut-off, we get terms that go as ^ as m tends to zero. This term multiplies the whole 
expression and does not drop out on differentiation. 



3. Going to de Sitter space Up to this point we have studied a metric for an impulsive 
wave propagating through minkowski space. We ran into problems which indicate that 
we may be perturbing around the wrong vacuum. Since the exact treatment of a related 
problem has no infrared divergences, we know that these divergences are spurious . 

As a possible way out we take the impulsive waves propagating in a de Sitter universe. 
Since the de Sitter space has a parameter with the dimensions of mass, one may think of 
trading thesec two parameters for one another. 



we try to perturb around the de Sitter space. Since the impulsive wave solution in 
the de Sitter space [11] is just the Minkowski solution multiplied by the factor (1 + Ap 
, we get the de Sitter two-point function just by multiplying the Minkowski case by (1 + 
A_p^ 4- Am 6 " ). This expression can be expanded as 

Auv\ ( Au'v'\ ( AUV\ 2 A, ... „ a2 , , 
1 + — ) ( 1+— J = + — 5— J + —(u-u')(v-v')+A 2 (...). (19) 
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Here U = ^ , V = ■ We multiply the Green function obtained in Minkowski space 
by this expression to get the de Sitter expression. In the Minkowski space is given above 
gives rise to Gp whose first term goes as 

®Ml + e (*V 2 (20 ) 

(s — s') 2 m 2 (s — s') 2 m 2 

for x=x',y=y'. Here m 2 is an infrared parameter. We introduced the infrared parameter 
by adding 2m 2 u 2 to our operator in equation (4). When the calculation is done in the 
presence of this additional term, we get the expressions given in equation (20) when we 
perform the summation operation to obtain the propagator function. In the presence of 
the new term the operators given in expression (18) are modified and result in terms given 
in (20) . When we go to de Sitter space we have to multiply them by the expression given 
in (19). 

Note that going to de Sitter space was only a technical trick. We will take A equal to 
zero at the end of the calculation and land in Minkowski space. We see that many of the 
terms that part of Minkowski space is subtracted for regularization are set to zero when 
A goes to zero. The terms given above are undeterminedsince they are multiplied by -^j, 
where both m 2 and A tend to zero. At this point we choose A proportional to m 2 . This 
choice is dimensionally correct. Any other choice, say, A proportional to the first power 
ofm times s, would be unnatural. To obtain the vacuum expectation value of the stress- 
energy tensor, we have to differentiate the Green function with respect to the coordinates. 
We are particularly interested in < T Mi , >. We first take m and A going to zero limit, and 
then differentiate with resoect to v and v' symmetrically, and then take the coincidence 
limit. We see that at the end of this calculation < T Mi , > turns out to be proportional to 

^4^. Here S is the Dirac delta function. The proportionality constant between the scalar 
curvature and the infrared mass squared can be absorbed in the perturbation constant a 2 . 



Conclusion We tried to calculate here the quantum fluctuations resulting from snapped 
cosmic strings, by perturbing around the vacuum, for a warp function that corresponds 
to the Gleiser-Pullin solution [14]. we could not seperate a finite part in the first-order 
calculation, in the second order, we ran into infrared divergences. If we use the same warp 
function in the calculation where we perturb around the de Sitter space and take the scalar 
curvature A proportional to m 2 , we can get finite results in second-order perturbation 
theory. In the above expression m is the infrared cut-off. At the end we get < > 



proportional to This is just the result we anticipated in the Minkowski case. At the 
end we are back in Minkowski space. 

This result suggests that to get finite results in perturbation theory for the case studied 
we have first to break the conformal invariance which does not allow vacuum fluctuations 
for the massless particles [10], by hand, by going to de Sitter space, and then come back to 
Minkowski space after the de Sitter space parameter cancels the infrared divergence. As a 
result of this operation we get a non-zero contribution for vacuum fluctuations. 

We anticipate somehow detecting the presence of cosmic strings [15]. We were ex- 
pecting to get a finite vacuum expectation value for the stress-energy. We could not get 
finite contributions in our previous calculations when we stayed in Minkowski space [4-7]. 
When similar calculations were performed for the de Sitter case, we got finite results which 
were proportional to positive powers of the curvature of the de Sitter space. [12,8]. The 
contributions were zero when the curvature was set to zero. 

In second-order calculations we encountered infrared divergences, which made the ex- 
pansion ambigious at this order. We cam give a meaning to these calculations by first 
starting with the de Sitter impulsive wave where the curvature of the space is proportional 
to the square of the mass of the scalar field whose stress-energy is calculated in the back- 
ground metric of the impulsive wave. This mass is used as an infrared parameter further 
in the calculation. 

Our method may be criticized because it was applied only to a single choice for the 
warp function, although our choice was a very important case, we tried to perform the 
same calculation for the case where 



HO = ) (21) 

and expanded the operator in powers of e. This choice is known to behave exactly the 
same way as the more general case, 

the second-order equation reads 

(L - A o )0 (o) = -L^°) - L'rfW (23) 

where 

T , 8v ( d d \ ^ 3d, 

L ^ = ^{ V Yv- U ^u) +F[v ^ X ^Yy ] - (24) 

Here O , 4>i an d 02 are the zeroth-, first- and second-order solutions. The first term on the 
right-hand side of this expansion gives terms that read as 
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which are of the same form as the terms in equation (16). When these terms are integrated, 
we obtain 

* ,2) - ^ {4 + m - ^ - + - (26 > 

where dots represent the additional terms we get from those included in F in equation 
(24). 

Doing the same calculation to obtain Green function, we get an expression that be- 
haves exactly as those given in equation (20). There are no terms which cancel these terms. 
At the end < T Mi , > comes out to be proportional to The details of this calculation 
will be published elsewhere. 

We think that this second calculation shows that our result is not restricted to one 
particular case, but illustrates the general behaviour in this problem. By starting with a 
massive field in the background metric of the impulsive wave in de Sitter space, we are 
breaking conformal invariance by hand. In the absence of the impulsive wave, i.e. in the 
zeroth-order perturbation theory, we do not have conformal invariance in the model, in 
contrast to the case when we have a massless scalar field and an impulsive wave in the 
minkowski space. Then we set the two parameters proportional to each other and send 
them to zero. As a result of this operation we get a finite contribution for < T^ v > . 

Our problem did not have conformal invariance in the first place. Using perturbative 
methods we could not convey this information to our solutions. We propose , then, first 
introducing additional parameters where this invariance is broken by hand in the pertur- 
bative calculation. At the end we send these additional parameters to zero in thwe same 
manner. We found that the peerturbative calculation gave the same qualitative answer, 
at least, that an exact solution would have given. 
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